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METHOD OF WATERMARKING A COLOR IMAGE BY 
VECTORIAL APPROACH, METHOD OF DETECTING A MARK 
OF WATERMARKING, CORRESPONDING DEVICES, IMAGE 
AND COMPUTER PROGRAMS 

5 FIELD OF THE INVENTION 

The field of the invention is that of the 
watermarking of still color images. More precisely, 
the invention concerns a watermarking technique which 
allows signatures to be incorporated into an image, 
10 taking into account the color dimension of the 
latter. 

BACKGROUND OF THE INVENTION 
The field of image tattooing, which is 
still called watermarking, is currently booming, and 
15 is the subject of much research, both in the fields 
of video sequences and still images. 

We can currently distinguish two main 
families of image watermarking, respectively 
corresponding to the watermarking of gray levels 
20 images and of color images. 

Indeed, the first years of research in the 
field of watermarking were mainly based on techniques 
for protecting gray levels images, and led to the use 
of three distinct fields of marking: the fields of 
25 space, frequency and multi-resolution. 

The watermarking models in the space field 
(where the mark is applied directly to the pixel 
values) have the advantage of only requiring a short 
calculation time. They are generally resistant to 
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geometrical attacks (for example rotation and 
changing of scale) . The methods proposed are 
histogram modification (Coltuc D. et al., "Image 
authentication by exact histogram specification" , 
5 workshop on multimedia signal processing, Cannes, 
France, October 2001) or again patchwork techniques 
(D- Gruhl, W. Bender, Moritomo, "Techniques for data 
hiding" in processing SPIE, volume 2420, page 40, 
February 1995) . 

10 As concerns the watermarking models in the 

frequency field, they have the advantage of being 
resistant to compression (such as JPEG for example) . 
The mark is applied to the coefficients resulting 
from a Fourier type transformation (as presented for 

15 example by V. Solachidis et I. Pitas, " Self -similar 
ring shaped watermark embedding in 2-D DFT domain", 
10 th European Signal Processing Conference 
EUSIPCO' 2000, Tampere, Finland, pages 1977-1980, 
September 2000) or Discret Cosinus (as presented for 

20 example by F. Alurki and R. Mersereau, "A robust 
digital watermark procedure for still images using 
DCT phase modulation" , 10 th European Signal Processing 
Conference EUSIPCO' 2000, Tampere, Finland, pages 
1961-1964, September 2000). 

25 Marking in the multi-resolution field 

offers several advantages. It is above all the field 
used in the most recent compression standards. It 
also allows the frequency band which will bear the 
mark to be chosen, thus permitting the risks of the 
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image being damaged by the application of the mark to 
be reduced (as illustrated by D. Kundur and D. 
Hatzinakos, "Digital watermarking using multi- 
resolution wavelet decomposition" , Proceedings of 
5 IEEE ICASSP '98, vol. 5, pages 2969 - 2972, Seattle, 
WA, USA, May 1998) . 

Apart from these watermarking techniques 
for gray levels images, a second main family of 
watermarking techniques for still images proposes 
10 taking into account the color dimension of the 
images . 

In this family, we can first of all 
distinguish a first sub-family of techniques 
consisting of adapting the gray levels methods to the 

15 three color components. Parameters are then used to 
control the force of the marking on each component in 
order to take into account the characteristics of the 
human viewing system. 

A second sub-family groups methods which 

20 are specific to the color images. It takes into 
account the human viewing system and uses the 
characteristics of the color representations. 

For example, one particular method, 
specific to color images, is that proposed by L. 

25 Akarun, N. Ozdilek, B. U. Oztekin, "A Novel Technique 
for Data Hiding in Color Paletted Images", 
Proceedings of the 10th European Signal Processing 
Conference, EUSIPCO'00, Tampere, Finland, pages 123- 
126, September 2000. 
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The first step of this method consists of 
quantifying the color space using the 'median-cut' 
algorithm. This consists of iteratively partitioning 
the color space with planes that are perpendicular to 
5 the color axes and passing through the median values 
of the data. 

There are two possible cases: 
all of the palette is used in the 
representation of the image, or 
10 - certain values of the palette are not used by 

the image. 

The first case is then no longer adapted to 
the marking algorithm. However the palette obtained 
includes colors that the human eye cannot 
15 distinguish. By using this property, certain colors 
of the palette can be liberated, so that they are not 
used in the image to be marked. 

The author proposes that two colors are 
indiscernible if A£<3 where 

20 AE = ^j(L { -L 2 ) 2 +(a x -a 2 ) 2 + (b { -b 2 ) 2 in the Lab space. A 
reminder is made that the Lab representation space is 
a perceptually uniform space. L represents the 
luminance and the a and b components are chromatic. 

Let M(i) be the mark, composed of colors 

25 that are not part of the palette. The author points 
out that a binary mark (made up of two colors) is 
more robust (the risk of detection error is thus 
reduced) . 

Let C be the most frequently used color of 
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color histogram, each absciss corresponding to a 
color of the palette.) 

The mark contains a lower number of 
elements than the number of color pixels C, i<h(C). 
5 Each of these elements is indiscernible from the 
color C. 

The marking consists of replacing the i th 
pixel of color C by M(i) . 

Another original method is that proposed by 
10 S. Battiato, D. Catalano, G. Gallo, R. Gennaro, in 
"Robust Watermarking for Images based on Color 
Manipulation" , Proceedings of the 3rd Workshop on 
Information hiding, LNCS 1768, pages 302-317, 
Dresden, 1999. According to this method, the mark is 
15 not created for the image beforehand, but it is the 
colorimetric content of the image which represents 
the mark. One disadvantage of this technique is that 
it requires a lot of data to be stored. 

The color space proposed by the author 
20 respects two properties: 

- the space must be perceptually uniform (as are 
the Lab and Luv spaces) so that a Euclidian distance 
measurement can be assimilated to a difference in 
color for the human viewing system; 
25 - the switch to this space, noted IC,C 2 , must be 

rapid, simple and have no loss of information. 

The field of opposed colors is defined here 
from RGB, as follows: 
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RGB^LC X C 2 : 
L = R + G + B; 

<C, =2B-R-G; 
C 2 =R-2G + B. 



LC X C 2 ^>RGB: 

R = (L + C 2 -C,)/3; 
<G = (L-C 2 )/3 

B = (C } +L)/3. 



According to the authors, this space is the 
closest to the representation of the chromatic 
5 channels of the human viewing system. 

The marking is carried out as follows. The 
mark is a vector M(n) , n=l,...,k,...,N where N is the 
number of colors of the image and Jc is a color of the 
image. Let (LC x C 2 ) k be the color associated to the 

10 index k, represented by a vector in the LC X C 2 space. A 
radius is selected in a sphere defined around the 
coordinate point (LC x C 2 ) k randomly. For each pixel 

(x,y) corresponding to the color k, the vector 
corresponding to the (previously determined) radius 
15 is added to the initial color vector to obtain the 
marked vector {LC x C 2 ) k . 

Each color is thus marked by adding a same 
vector. The marked image is reconstructed by 
replacing the original colors with the marked colors, 
20 by respecting the coordinates of the image pixels. 

The mark is then formed by all of the 
colors of the original image. 

Detection is made by comparing the marked image with 
the marks generated on all of the images processed. 
25 This is carried out as follows. Firstly, it must be 
considered that the image on which the detection is 
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to be made has the same number of colors as the 
original image. The detection algorithm compares the 
image to each of the recorded marks. Color by color, 
a search is made for the closest mark (where the 
5 difference vectors between the two pixels of the same 
coordinate are the most similar) . The mark 
corresponding to the marked image is that with the 
highest number of difference vectors between the 
colors of the marked image and the original image in 
10 common. 

It can be noted that the various 
publications concerning color watermarking techniques 
are generally consecrated to the integrity of the 
visual appearance and the robustness of these 
15 techniques against classic attacks is hardly 
developed. 

Generally, for methods whose basic 
algorithm can be applied to gray levels images, the 
marking uses the following technique: 
20 - switch in the transformed space (wavelet 

coefficients, discrete cosine coefficients, etc.) 
- application of the following formula: 

I w (iJ) = r(iJ) + a{iJ)M{iJ) 

where I w represents the transform of the marked image 

25 (or component), I' the transform of the original 
image (or component) , M the mark, a the marking 
intensity control factor, and where i and j represent 
the coordinates of the processed pixel. 

The algorithms based on the characteristics 



-8- 



of the color components are more adapted to the 
characteristics of the human viewing system. The 
first article of Kutter M., Jordan F. and Bossen F, 
{"Digital Signature of Color Images using Amplitude 
5 Modulation", Processing of SPIE storage and retrieval 
for image and video databases, San Jose, USA, volume 
3022, number 5, pages 518-526, February 1997) on 
color watermarking proposed working on the blue 
component of the RGB system, for which the human eye 

10 is the least sensitive. 

This idea was also taken up by A. Reed and 
B. Hannigan, in "Adaptive Color Watermarking" , 
Proceedings of SPIE, Electronic Imaging, volume 4675, 
January 2002. These authors propose working on the 

15 yellow component of the system CMY (Cyan Magenta 
Yellow, which is a colorimetric space) , as they 
considered that the human eye is less sensitive to 
color variations on the yellow - blue axis. 

Generally speaking, the blue component 

20 allows the mark to be hidden more effectively (but it 
is not as strong) and the green component allows the 
mark to be protected more effectively against attacks 
(but with greater visual damage to the image) . This 
is explained by the fact that the human viewing 

25 system is more sensitive to variations in green than 
in blue. The compromise between invisibility and 
strength of the mark therefore depends on the 
colorimetric nature of the component. 

Finally, a last technique, proposed by J.J. 
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Chae, D. Mukherjee, and B.S. Manjunath, in "Color 
Image Embedding using Multidimensional Lattice 
Structures" , Proceeding of IEEE International 
Conference on Image Processing, Chicago, Illinois, 
5 volume 1 , pages 460-464, October 1998, is based on a 
vectorial type approach. 

According to this technique, the first step 
of the marking consists of breaking down into 
wavelets the initial image and the signature (which 
10 may also be an image) . A single breakdown level is 
used. In this way, the wavelet coefficients of the 
original image are obtained, noted (C Y , C v , C v y) , for 

which each component corresponds to the color 
component of the YUV space, and the wavelet 
15 coefficients of the mark. The use of the YUV space 
(space used for video: Y is the luminance component, 
U and V are chromatic components) thus allows a 
direct adaptation of this technique to video 
documents . 

20 The wavelet coefficients resulting from 

this breakdown of the mark are quantified in /? 

levels. We thus obtain a vector M , (M Y , M u , M v y) 

for a color mark and M(x,y) for a gray levels, 

containing (j.) elements where 1 </</?. 

25 The integration of the mark can be written 

as follows: 

(C Y X u Xy)i\x,y) = (Cy,C u ,C y ) i (x 9 y) + aM(s i ) 
where a is the control factor for the marking force. 
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The detection of the mark is then carried 
out as follows. After application of the transform in 
wavelets onto the marked image, the resulting 
coefficients are quantified in levels. 
5 To estimate the closest vector to that used 

for the mark, a search is carried out for the 
quantification value of the closest coefficient to 
those of the initial image (the mark is then 
detected, element by element) . 
10 The image corresponding to the mark 

detected can thus be reconstructed by transformation 
of the inverse wavelets. 

All of the gray levels or color 
watermarking techniques described above have a number 
15 of disadvantages. 

Therefore, the watermarking techniques for 
color images based on gray levels algorithms do not 
permit the color dimension of the images to be taken 
into account. Consequently, they do not permit the 
20 compromise between the invisibility and the strength 
of the mark to be optimized. 

The color image watermarking techniques are 
very complex and all attempt to satisfy an 
invisibility objective for the mark. In other words, 
25 they are not strong enough to resist classic attacks 
such as compression or filtering. 

The technique proposed by Chae et al. 
described above, even though it uses a vectorial 
approach, has the disadvantage of not being resistant 
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to compression. Among others, this marking method 
does not allow making documents secure to be 
envisaged. Finally, this technique does not take into 
account the color dimension of the image, as it 
5 attributes the same marking force to the three color 
components of the image. 

Finally, none of the techniques of the 
prior art using a wavelet breakdown offers blind 
marking, permitting detection of the mark without the 
10 original image being required. Furthermore, none of 
these techniques offers resistance to JPEG 
compression, median filtering and the addition of 
noise . 

The particular aim of the invention is to 
15 overcome these disadvantages of the prior art. 

More precisely, one aim of the invention is 
to provide a color watermarking technique permitting 
one or more signatures to be incorporated to a color 
image . 

20 In other terms, one aim of the invention is 

to provide a technique for making documents secure by 
the invisible insertion of a mark. 

Another aim of the invention is to 
implement such a technique that is resistant to most 

25 attacks. In particular, one aim of the invention is 
to provide such a technique which is resistant to 
compression (especially JPEG type) , to filtering 
(especially by a median filter) and the addition of 
noise . 
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The invention has another aim of 
implementing such a technique which permits blind 
detection, which is to say that does not require the 
original image, of the mark in an image. 
5 Another aim of the invention is to provide 

such a technique permitting the visibility of the 
mark to be reduced in comparison to the techniques of 
the prior art. 

These aims, as well as others which will 
10 become clearer further on, are achieved by a color 
watermarking process with at least three components. 

SUMMARY OF THE INVENTION 
According to the invention, such a process 
comprises a step where the watermark is inserted, on 
15 at least one point of the image, according to an 
insertion rule taking into account the relative 
position of at least three component vectors, for 
each of said components respectively, associated to 
said at least one point. 
20 In this way, the invention is based on a 

completely new and inventive approach of color images 
watermarking. In fact, the invention is based on a 
vectorial approach, permitting, in comparison with 
the techniques of the prior art, to take into account 
25 the color dimension of the image. 

The watermark is therefore no longer, 
contrary to the techniques of the prior art, applied 
identically to the three color components of the 
image; in other terms, whereas, according to the 
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prior art, independent marking of the three component 
vectors of the image was carried out, according to 
the invention, these three vectors are now taken into 
account in synergy for the insertion of the mark. 
5 The implementation of a vectorial approach 

permits, in terms of resistance, to obtain more 
stable responses than with the techniques of the 
prior art in the field of wavelets, such as that of 
Kundur et al. mentioned previously, as the 

10 watermarking technique in the wavelet field. Such a 
vectorial approach also permits the variations 
induced by the mark in the color field to be 
minimized, and therefore slightly improve the 
invisibility of the mark with respect to the 

15 techniques of the prior art. 

Advantageously, for each of the said points 
concerned, two vectors are selected as reference and 
one vector to be marked in order to bear the said 
watermark. 

20 It can be noted that the invention consists 

therefore, contrary to the prior art, of marking the 
component vectors (i.e. color vectors for example), 
and not just coefficients of the image. 

Preferably, the said vector to be marked is 

25 the central vector, situated between the said 
reference vectors. 

Advantageously, the distances are 

calculated two by two between said three vectors, 
said reference vectors being the furthest apart from 
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one another. 

Preferably, a frontier is determined 
between the said reference vectors, defining two 
zones respectively associated to the binary values 
5 "0" and "1". 

Two half spaces are thus created, 
respectively associated to each of said reference 
vectors, and corresponding to a value "0" or "1" of 
the mark that is to be inserted. 
10 Advantageously, said frontier is defined by 

the bisector between said reference vectors. 

According to one advantageous 

characteristic of the invention, the marking of the 
said vector to be marked comprises a possible shift 
15 of the said vector to be marked to one of the said 
zones, according to the binary value of the marking 
to be applied. 

Therefore, if a value "0" mark is to be 
inserted on the vector to be marked and the vector to 
20 be marked is in the zone or the half space associated 
to the binary value "1", its end is moved so that it 
is positioned in the other half space, or zone, 
associated to the value "0". 

If a value "0" mark is to be inserted on 
25 the vector to be marked and the vector to be marked 
is already in the zone associated to the binary value 
"0", its end may also be moved so that it is moved 
closer to the reference vector located in this half 
space, in order to increase the resistance of the 
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mark (while remaining of course in the half space 
associated to the binary value "0") . 

According to another advantageous 
characteristic of the invention, the said movement 
5 has a variable amplitude, according to a marking 
force selected to suit at least one local 
characteristic of the said image. 

In this way, according to the colorimetric 
characteristics or local texture characteristics of 

10 the image, a mark that is more or less visible and 
more or less resistant can be chosen, according to 
the marking force selected. Such a marking force may 
adopt for example a value included between 0 and 1, a 
nil force corresponding to a mark with low resistance 

15 and invisible, and a force equal to 1 corresponding 
to a strong but visible mark. 

According to one advantageous variant of 
the invention, such a watermarking process also 
comprises a wavelet transformation step for each of 

20 the said components of the image, and the said three 
component vectors are determined, for each point of 
at least one level of breakdown of the said wavelet 
transformation, for each of the said components 
respectively. 

25 The invention therefore combines with the 

vectorial approach a breakdown into wavelets of the 
image. The use of such a wavelet transformation 
advantageously permits the technique of the invention 
to be more resistant to JPEG type compression 
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attacks, or median filtering or the addition of 
noise . 

It can be noted that by "point it is 
meant here a pixel position in a level of breakdown 
5 of the image, as illustrated by figure 10, which 
represents an image scale E 101, and the three 
detailed images D N 1,R 102, D N 2 ' R 103, D N 3 ' R 104 for the 
red component of a wavelet breakdown in an RGB 
system. The 105 " point " is defined by a same pixel 
10 position in each of the three detailed images 
referenced 102 to 104. 

Preferably, at least two vector marking 
conventions are provided. 

According to an advantageous characteristic 
15 of the invention, the marking convention selected for 
a given image is that limiting the risks of conflicts 
when detecting the said mark. 

Preferably, one of the said conventions is 
selected according to the number of reference vectors 
20 present in the said image for each of the said 
components . 

Thus, in an RGB type system for example, 
the number of times that the red, green and blue 
vectors are chosen as reference or marker vectors is 
25 determined, and the marking convention permitting the 
risks of conflict during detection of the mark to be 
minimized is deduced. For example, a marking 
convention is selected so that the vector the least 
often marked is situated in the conflict zone. 
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Preferably, the said marking is, at least 
in normal situations, calculated using the following 
equation: 

v M A*>y) = fay) - 0 - f m W*(*,y) - K <* jO) 

where V M {x,y) is the said vector to be marked, V R (x,y) 
is one of the said reference vectors, F M is the said 
marking force and V MW (x,y) is the said marked vector. 

Preferably, in the event of a conflict, the 
said marking is calculated using the following 
equation: 

where jbeta<2. 

Advantageously, such a watermarking process 
comprises, after the said step of associating a 
watermark, and an inverted wavelet transformation 
step, providing a marked image. 

In this way, after carrying out a wavelet 
transformation of the various color components of the 
image, and after having inserted the watermark on the 
color vectors, a marked image is reconstructed by 
inverted wavelet transformation. 

Advantageously, the said mark is a pseudo- 
random binary signature written redundantly. 

The redundancy of the insertion of the mark 
permits better protection of the image and a higher 
probability that the watermark will be detected. 

Preferably, the said components belong to 
the group comprising: 

- the RGB components; 
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- the YUV components; 

- the CMY components. 

More generally, all of the spaces where the 
color is represented can be sued, whether it be the 
5 spaces defined by the CIE (International Lighting 
Commission) or original spaces. 

The invention also concerns a device for 
watermarking a color image, with three components, 
and comprising means for inserting a watermark, on at 

10 least one point of the image, according to an 
insertion rule taking into account the relative 
position of three component vectors, for each of the 
said components respectively, associated to the at 
least said one point. 

15 The invention also concerns a computer 

programme comprising programme code instructions 
recorded on a support that may be used in a computer 
for watermarking a color image, with three 
components. According to the invention, the said 

20 programme comprises programming means that may be 
read by the computer to carry out an insertion step 
for a watermark, on at least one point of the image, 
according to an insertion rule taking into account 
the relative position of three component vectors, for 

25 each of the said components respectively, associated 
to the said at least one point. 

The invention also concerns a color image 
with three components and comprising a watermark 
obtained by insertion of a watermark, on at least one 
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point of the image, according to an insertion rule 
taking into account the relative position of three 
component vectors, for each of the said components 
respectively, associated to the said at least one 
5 point. 

The invention also concerns a watermark 
detection process for a watermarked image, made 
according to the previously described watermarking 
process. According to the invention, such a detection 

10 process comprises a step for recovering a watermark, 
on at least one point of the image, according to a 
recovery rule taking into account the relative 
position of three component vectors, for each of the 
said components respectively, associated to the said 

15 at least one point. 

Preferably, the said recovery step 
comprises a sub-step for calculating the distances 
two by two between the said three vectors, and the 
said two vectors the furthest away from one another 

20 are the reference vectors, the third vector is the 
marked vector bearing the said watermark. 

Advantageously, a frontier between the said 
reference vectors defining two zones respectively 
associated to the binary values "0" and "1", the said 

25 recovery step also comprises an identification sub- 
step for the zone in which the said marked vector is 
situated, and a step for associating a corresponding 
binary value. 

Advantageously, the said mark is recovered 
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at least twice, and a correlation calculation is used 
with respect to a reference signature, in order to 
decide if the watermark has been detected correctly 
or not. 

5 The invention also concerns a device for 

detecting a watermark in a watermarked image, made 
according to the watermarking process previously 
described. According to the invention, such a device 
comprises means for recovering a watermark, on at 

10 least one point of the image, according to a recovery 
rule taking into account the relative position of 
three component vectors, for each of the said 
components respectively, associated to the said at 
least one point. 

15 The invention also concerns a computer 

programme comprising programme code instructions 
stored on a support that may be used on a computer to 
detect a watermark in a watermarked image, made 
according to the process previously described. Such a 

20 programme comprises programming means that may be 
read by a computer to carry out a recovery step for a 
watermark, on at least one point of the image, 
according to a recovery rule taking into account the 
relative position of three component vectors, for 

25 each of the said components respectively, associated 
to the said at least one point. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Other characteristics and advantages of the 
invention will become clearer upon reading the 
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following description of a preferred embodiment, 
provided simply as a non-restrictive illustration, 
and the appended diagrams among which: 

figure 1 shows a block diagram of the 
5 multi-resolution analysis principle for an image I by 
wavelet transformation, implemented according to the 
invention; 

figure 2 illustrates the principle for 
obtaining images at resolution j level from the image 
10 scale at level j+1 by wavelet transformation, 
according to the principle of figure 1; 

figure 3 describes a structure of wavelet 
coefficients formed after the wavelet transformation 
of figure 2; 1 shows a block diagram of the multi- 
15 resolution analysis 

figure 4 illustrates the generation of a 
mark by redundancy of the signature; 

figure 5 presents the construction step of 
the reference vectors and the marking vector; 
20 figures 6a to 6c illustrate the various 

marking conventions used according to the inventions- 
figure 7 illustrates one possible example 
of the position of the original vectors, as well as 
the two marking vector movements that can be 
25 envisaged for this case; 

figures 8a and 8b respectively present an 
example of an original image and a corresponding 
image marked using the watermarking technique of the 
invention; 
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figures 9a to 9c present, in the form of 
graphs, the results of resistance tests of the 
watermarking technique of the invention, on the 
example of the image of figure 8b; 
5 figure 10 states the notion of "point" of 

an image, which represents, as part of a wavelet 
transformation, a pixel position in a breakdown 
level . 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
10 The general principle of the invention is 

based on taking into account the color dimension of 
the images, and proposes a vectorial approach of 
watermarking. The invention thus permits vectorial 
marking taking into account, in synergy, the three 
15 color vectors of the image. 

In one specific embodiment of the 
invention, which will be described in all of the rest 
of the document, this vectorial approach is combined 
with a wavelet breakdown of the color components of 
20 the image. 

In this specific embodiment, the invention 
therefore consists of inserting a watermark into a 
color image by carrying out the following steps: 

- Transformation of each component of the color 
25 image in the wavelet field, 

- Pseudo-random independent generation of the 
mark based on a key, 

- Definition of all of the color vectors 

- For each coordinate of the passing band 
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components : 

Definition of 2 reference vectors and 
one vector which will bear the watermark 
(the reference vectors are the vectors 
5 which are the furthest away, the mark is 

thus applied to the "central" vector, which 
minimizes the visual impact) 

Definition of a rule which controls 
the marking force and determines the 
10 insertion space associated to the bit 1 and 

bit 0. 

- Reconstruction of the color image marked by 
the wavelet transformation, taking into account the 
modifications by marking of the wavelet coefficients. 

15 These various operations are described in 

greater detail further on in the document, in 
relation to the figures. 
1. Wavelet transformation 

The first step implemented according to the 

20 invention is a wavelet transformation step. A 
reminder is made that the wavelet transformation is a 
powerful mathematical tool allowing a multi- 
resolution analysis of a function, as described by 
Mallat S. (in "A Theory for Multi-resolution Signal 

25 Decomposition: the Wavelet Representation" , IEEE 
Transactions on Pattern Analysis and Machine 
Intelligence, Vol. 11, No. 7, July 1989, pp. 674- 
693), by Stollnitz E.J., DeRose T.D., et Salesin D. 
(in "Wavelets for Computer Graphics: A Primer-Part 
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2", IEEE Computer Graphics and Applications, May 
1995, pp. 76-84 or again in "Wavelets for Computer 
Graphics: A Primer-Part 2", IEEE Computer Graphics 
and Applications, July 1995, pp. 75-85). 
5 Appendix 1, which is an integral part of 

this patent application, makes a reminder of the 
essential aspects of the wavelet theory. 

In the scope of the invention, the 
functions considered are digital images, which is to 
10 say discrete bi-dimensional functions. With no loss 
of general features, it is supposed that the 
processed images are sampled on a discrete grid of n 
lines and m columns and at a value in a luminance 
space sampled at 256 values. Furthermore, it is 

15 supposed that n = 2 k (keZ) and that m = 2 ! (leZ). 

If we note that I the original image, this 
therefore gives: 

j0,m]x[0,n]-> [0,255] 
' (x,y)h+I(x,y) 

The wavelet transformation of the image I 
20 10 makes possible a multi-resolution representation 
of I, as illustrated by figure 1. At each level of 
resolution 2 J (j<-l), the representation of I 10 is 
provided by a rough image A^I 11 and by three 

detailed images D[jl 12, D\j 13 and D\jl 14. Each of 

25 these images is of the size 2 k+J x2 l+J . 

Again, at the following level of 
resolution, the rough image A J 11 may be represented 
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by a rough image A^.J llx and by three detailed images 
& 2 „I 12 x , D\ } J 13i and D\,J 14 x . 

The wavelet transformation requires the 
selection of a scale function <&(x) as well as the 
selection of a wavelet function ^(x) . Based on these 
two functions, we derive a scale filter H and a 
wavelet filter G whose respective pulse responses h 
and g are defined by: 

h(n) = (# rl (u),#(u-nj)VneZ 
g( n) = (i// vl ( u ),</)( u - nfpfn e Z. 

Let us note respectively H and G the 
mirror filters of H and G (i.e. h(n) = h(-n) and 

g(n) = g(~n)) . 

We can therefore show, as illustrated by 
figure 2, that: 

- A 2 jl 20 can be calculated by convoluting A 2/+ J 

21 with H in both dimensions and by sub-sampling a 
factor two (23, 24) in both dimensions; 

- D ] 2J I 22 can be calculated by: 

1. convoluting A 2j+l I 21 with H in the y 

direction and by sub-sampling a factor two (23) in 
this same direction; 

2. convoluting the result of step 1) with G 
in the x direction and by sub-sampling a factor two 
(25) in this same direction. 
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- D^jl 26 can be calculated by: 

1. convoluting A^J 21 with Gin the y 

direction y and by sub-sampling a factor two (27) in 
this same direction; 

5 2. convoluting the result of step 1) with H 

in the x direction and by sub-sampling a factor two 

(28) in this same direction. 

- D 3 2J I 30 can be calculated by: 

1. convoluting A^J 21 with G in the y 

10 direction and by sub-sampling a factor two (27) in 
this same direction; 

2. convoluting the result of step 1) with G 
in the x direction and by sub-sampling a factor two 

(29) in this same direction. 

15 In one particular embodiment of the 

invention, the Daubechies base is used, and the level 
of resolution 2 r (r<-\) where r=-4 is selected. 

The Daubechies wavelets have a minimum 
sized support for a given number of nil moments. The 

20 compact support Daubechies wavelets are calculated 
from conjugated mirror filters h with finite pulse 
responses. The length of the filter determines the 
time support of the associated wavelet as well as the 
number of nil moments. There is no explicit formula 

25 for the Daubechies filters regardless of the order of 
the filter. The filters are constructed based on the 
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resolution of the following problem: 



10 



AM 



Let h{w) = Yj h \P\ e ' inW be the Fourier 



n=0 



transformed trigonometric polynomial of the 
Daubechies filter, therefore a polynomial R{e^) has 
to be constructed with a minimum degree such that 



h(w) = V2 



Rie'™) and h(w) 2 + h(w + n)\ 2 =2. According 



2 

to the invention, the filter with eight coefficients 
is preferred, which corresponds to the following 
polynomial : 

h(w) = 0,2304 + 0,7148*-*" + 0,6309.<f ' 2w -0,028if -0,187iT'" 4w 
+ 0 5 0308^ ,5w + 0 3 0329if ' 6w - 0,01 06£ nw 
Figure 3 illustrates the construction, 
after the wavelet transformation has been carried out 
to a level of minimum resolution 2 r (r<-l), of a 
wavelet coefficient structure: 
15 - Each pixel p(x,y) referenced 34 on the image 

A y I is the root of a tree structure; 

- Each root p(x,y) is attributed three thread 
knots referenced 31 to 33 designated by the wavelet 
coefficients on the three detailed images D S 2 J 

20 (s=l,2,3) located in the same position (x,y) ; 

- Due to the sub-sampling using a factor two 
carried out by the wavelet transformation at each 
change of resolution (see figure 2), each wavelet 

coefficient D s r (x,y) (s=l,2,3) corresponds to a zone 
25 whose size is 2x2 pixels in the detailed image 
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corresponding to the resolution 2 r+ . This zone is 
located in (2x,2y) and all of the wavelet coefficients 

belonging to it become thread knots of D s r (x,y). 

Recursively, the tree structure is thus 
5 constructed in which each wavelet coefficient D* 9 (x 9 y) 
(5=1,2,3 et 0>u>r) has four thread knots designated 
by the wavelet coefficients of the image D s r +J located 
in the region situated in (2x,2y) whose size is 2x2 
pixels . 

10 Once the tree structure has been 

constructed, each wavelet coefficient D s r (x,y) 

(s=l,2,3) corresponds to a region whose size is 2~ r x2~ r 

pixels in the detailed image D* r J . 

2. Construction of the mark 
15 A binary S signature of N 2 bits is 

generated in a pseudo-random manner, checked by a key 
K. This signature, written in the form of a matrix 

Vtf is then written in a redundant manner to 

obtain the mark W. The redundancy can be carried out 

20 bit by bit or signature by signature, as is 
illustrated in figure 4. 

In a particular embodiment of the 
invention, the signature S 40 is a 16 bit binary 
signature and is generated in a pseudo-random manner 

25 using a key K=9. This bitmap signature (4*4) 40 is 
then copied T times to create a mark, as in the first 
example 41 shown in figure 4. For an image of size 
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256*256, knowing that the wavelet coefficients 
employed are those defined for the 4 th scale of 
decomposition, the mark will be of size 16*16. The 
signature is therefore repeated 16 times. 
5 We note that in figure 4, the black squares 

illustrate a bit by 1 and the white squares a bit by 
0. 

3. Insertion of the mark 
3.1 Definition of the vectors 
10 The vectors are defined from the 

decomposition into wavelets at the scale -r, for each 

co-ordinate (x,y) of the three detailed images D^l 

22, D\ } I 26, D\jl 30, and for each of the three 

designated colored components I={A, B, C}. 
15 A, B and C are components that depend on 

the representation space (ABC can be RGB, XYZ, YUV, 
etc. ) . 

We thus obtain, for each co-ordinate (x,y) 
50 of the detailed images, the following three 
20 vectors, illustrated in figure 5: 

V A (x, y) = {D 2i A(x, y) ? D*j A(x, y), D 3 2j A(x, y) } , 
V B (x, y) = {D l 2j B(x, y), D^x, y), D 3 2i B(x, y) } , 
V c (x, y) = {D 2j C(x, y), D^C(x, y), D 3 2j C(x 5 y) } . 

In a particular embodiment of the 
25 invention, we examine the image in the representation 
space RGB, and the vectors are defined for each 
component (R, G and B) from the wavelet coefficients 
of the 4 th scale of decomposition (r=4) : 
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V R (x,y) = {D'^Rfcy), D^R(x,y), D^R(x,y) } , 
V G (x 5 y) = { D;. 4 G(x, y), G(x, y), D 3 ^ G(x, y) } , 
V B (x, y) = {D;. 4 B(x, y), B(x, y), D 3 2 , B(x, y) } 
3.2 Calculation of distances, reference vectors and 
marked vectors 

For each co-ordinate of the detailed 
images , the distances between any two vectors are 
calculated: 

D A . B (x,y) = |v A (x,y)-V B (x 5 y)| 
D A; c(x ) y) = |v A (x 5 y)-V c (x J y)| 

DB, c(x 5 y) = |v B (x,y)-V c (x 5 y)| 

The greatest distance is defined by the two 
reference vectors V Rj 51 and V Rz 52 (see figure 5) . It 

is the third vector V M 53 which will be marked. 

Thus, if D AB (x,y)>D AC (x,y) and 

D^ B (x,y)>D B)C (x,y), then V A (x,y)and V B (x,y) are the 

reference vectors and V c (x ? y) is the vector which will 
bear the mark. 

If D AC (x, y) > D A B (x, y) and D A c (x, y) > D BC (x, y) , 

then V A (x,y) and V c (x,y) are the reference vectors and 

V B (x,y) is the vector which will bear the mark. 

If D BC (x,y) >D AB (x,y) and D BC (x,y) >D AC (x,y) , 

then V B (x,y) and V c (x,y) are the reference vectors and 

V A (x,y) is the vector which will bear the mark. 

In a particular embodiment of the 
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invention, the distances between the vectors are 
calculated, between any two, for each of the 256 
possible co-ordinates (16*16 wavelet coefficients in 
the detailed images) . 

D RiG (x,y) = |v R (x,y)-V G (x,y)| 
DR.B(x,y) = |v R (x,y)-V B (x,y)| 

DG,B( x 'y) = |^G(^y)-v B (x jy )| 

If D RG (x,y)>D RB (x,y) and D R G (x,y) >D G B (x,y) , 

then V R (x,y) and V G (x,y) are the reference vectors and 

V B (x,y) the vector which will bear the mark. 

If D RB (x,y) >D RG (x,y) and D R B (x,y) >D G B (x,y) , 

then V R (x,y) and V B (x,y) are the reference vectors and 
^c( x '^) vector which will bear the mark. 

If D GB (x,y)>D RG (x,y) and D GB (x,y) >D RB (x,y) , 

then V G (x,y) and V B (x 9 y) are the reference vectors and 

V R (x,y) the vector which will bear the mark. 

3.3 Definition and choice of the marking agreement 

For each co-ordinate of the detailed 
images, we try to know the number of times the 

vectors V A (x,y) , V B (x,y) and V c (x ? y) are used as 

references. This stage thus allows to choose the best 
agreement for the insertion operation (see figures 6a 
to 6c) , in order to minimize the conflict which will 
be explained in fuller detail in §3.5. 

If V A (x,y) is the least used vector as a 
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reference, the choice will be for agreement 
3, illustrated in figure 6c. 

If V B (x 5 y) is the least used vector as a 
reference, the choice will be for agreement 1, 
5 illustrated in 6a. 

If V c (x 5 y) is the least used vector as a 

reference, the choice will be for agreement 2, 
illustrated in figure 6b. 

Thus, in case of the "Home" image 
10 illustrated in figures 8a (original image) and 8b 

(marked image) , the vector V R (x,y) is used 235 times as 

a reference vector, the vector V G {x,y) 108 times and 

the vector V B {x,y) 169 times. To minimize the 
conflicts, we therefore use agreement 1, shown in 
15 figure 6a. 

However, it is recommended to make a 
precise choice in the case of small images. The 
automatic choice of the marking agreement is better 
adapted to large size images (at least several 

20 thousand pixels) . 

3.4 Insertion operation 

Marking consists in modifying the vector by 
placing its end nearer to the end of one of the two 
reference vectors, on the basis of the value of the 

25 mark (0 or 1), and this being for each co-ordinate 

50, as is illustrated in figure 5. We note V R] (x 9 y) and 
V R2 (x,y) (51 and 52) as the reference vectors and 
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V M (x,y) as the vector to be marked 53. 

Generally speaking, the marking is in the 

form of : 

5 with depending on the case: Vr=Vr\ or V R =V R2 
according to the value of M and the chosen agreement. 

F M represents the marking force applied to 
the vector with M={A, B, C}. The marking force can 
thus be adapted according to the local 
10 characteristics (colorimetric or texture) of the 
image . 

Two possible situations could develop, as 
is illustrated for a particular case in figure 7. 

Depending on the chosen agreement and the 
15 value of the mark W(x,y): 

- the vector V M (x,y) 73 is already in the correct 
half space where V MW (x 9 y) must be positioned. In this 
case, the formula (1) can be directly applied to 
V M (x 9 y) 73. This case is illustrated in figure 7, in 

20 consideration of the particular case where W(x,y)=l, 
and where the agreement imposes that V MW (x,y) is in 

the half space containing V R2 (x 9 y) 72. 

- the vector V M (x 9 y) 73 is not in the correct 
half space where V MW (x 9 y) must be positioned. In this 

25 case, prior to applying the formula (1) to V M (x 9 y) 73, 
the position of the end of this vector must be 
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modified. To do so, this end must be displaced in a 
collinear direction to the vector V coi (x 9 y)=V RX {x>y)- 
V R2 (x,y) , for a distance (D m +s) , where D m is equivalent 

to the distance between V M (x,y)13 and (P M ) (where P M is 
5 the border between the two half spaces associated to 
each of the reference vectors (71 and 72) . Preferably 

we choose e equals 10% of D m , so that V M (x,y) 73 does 
not get confused with (P M ) . This case is illustrated 
in figure 7, in consideration of the particular case 
10 where W (x,y) =0, and where the agreement imposes that 

V MW (x,y) is in the half space containing V Rl (x,y) 71. 

We note that the case where F M =0, the 
marking is of minimal robustness. The vector V M (x,y) 
73 becomes V M ^{x 9 y) 9 by simply displacing it (if 

15 necessary) in the half space corresponding to the 
positions of the reference vectors 71 and 72, the 
value of the mark and the chosen agreement. The image 
is thus degraded to the minimum. 

In the case where F M =\, the marking force 

20 is at the maximum. The mark is very robust, but has a 

high probability of being visible. V MW (x,y) thus 

becomes confused with one of the reference vectors 71 
and 72 on the basis of the previously indicated 
situation . 

25 The fact of choosing an intermediary value 

of F M allows to control the compromise between the 
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invisibility and the robustness of the mark. 

3.5 Conflict 

The marking agreement diagrams define the 

direction of the marking for each of the possible 
5 situations. In the case where the marking force F M =l, 

that being at the maximum, the marked vector becomes 

equal to a reference vector. 

Thus, in figures 6a to 6c, we note that, 

for each of the agreements, the bits of the mark are 
10 opposite for one of the branches of the agreement 

triangles. A conflict is thus possible. Indeed, 

during the transmission of the marked image or the 

processes it can undergo, the two confused vectors 

can be substantially modified. It is therefore 
15 possible that the reference vector becomes the marked 

vector and vice versa. 

During the detection of the mark (which 

will be described in fuller detail in §5 below) , 

detection errors can thus be integrated due to this 
20 conflict. The conflict is minimized if the chosen 

agreement minimizes this situation. The agreement 

choice stage proposed in §4.3 above is therefore 

particularly important . 

In the case of a conflict, the marking 
25 operation becomes: 

where beta<l. 

Thus the choice of the agreement allows to 
minimize the number of conflicts, and the errors due 
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to the conflicts are thus also minimized through 
applying a beta factor in the marking operation. 

In a particular embodiment of the 
invention, the marking forces employed are the 
5 following: F R =0 . 4 (for the red component); F G =0.4 (for 
the green component); F B =0.8 (for the blue component), 
with jbeta = 0.9 (in the case of a conflict). 
3.6 Reconstruction of the image 

The image is then reconstructed by 

10 transformation into inverse wavelets, taking into 
account the coefficients modified by the marking. 
Thus, figure 8b illustrates the marked "home" image, 
reconstructed by transforming into inverse wavelets, 
and figure 8a shows the same image prior to marking. 

15 4 . Detection of the mark 

From now on we will describe the different 
operations implemented when we try to detect a mark 
in a watermarked color image according to the 
invention. Indeed, after processing the image 

20 (transmission, JPEG compression, median filtration or 
noise addition) the mark can be detected. 

4.1 Decomposition into wavelets 

The decomposition operation in the field of 
wavelets of the marked image is identical to that 
25 carried out in the context of marking, shown in §1 
above . 

4.2 Definition of vectors and choice of agreement 

Likewise, the operations for defining 
vectors and the choice of marking agreement are 
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identical to those previously described in §3.1, 3.2 
and 3.3. 

4.3 Recuperation of a bit of the detected mark 

For each co-ordinate, we check the 
5 proximity between the marked vector and the two 
reference vectors. Depending on the chosen agreement, 
it is possible to recognize the value of the 
integrated bit. The table below shows all the 
possible situations in the reconstruction of the 
10 mark. The mark is thus reconstructed bit by bit. 





convention 1 


convention 2 


convention 3 


V A [x,y) - V B {x,y)\ < \v A (x,y) - V c {x,y)\ 


0 


0 


1 


v A {*,y) - v c (*,y)\ < \vAx,y)-v B {x,y)\ 


1 


1 


0 


V B (x>y) - VA(.*,y)\ < f B ^,y) - V c {x,y)\ 


0 


1 


1 


V B (*> y) - V c (x, y)\ < \V B (x, y) - V A (x, y)\ 


1 


0 


0 


V c (x,y) - V A (x,y)\ < f c {x,y) - V B (x,y)\ 


0 


1 


0 


y c (x,y) - V B {x,y)\ < \v c (x,y) - V A (x,y)\ 


1 


0 


1 



4.4 Reconstruction of the signature 



The detected signature is then 
reconstructed from the detected mark. To do so, the 
mean of the bits of the mark corresponding to a co- 
15 ordinate of the signature is carried out, according 
to a principle similar to that used during the 
construction of the mark described in §2. Thus the 

detected signature S is reconstructed, bit by bit. 

In other terms, we start, for each co- 
20 ordinate of the detailed images, by referring to the 
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table above, which indicates the value of the 
inserted bit, for each of the agreements in figures 
6a to 6c, according to the distances between the 
vectors . 

5 In a particular embodiment of the 

invention, the signature S of 16 bits is the 
reconstructed by calculating the mean of the 16 
values of the mark for each of the corresponding co- 
ordinates, according to the form of the redundancy 

10 defined in §2. 

4.5 Recognizing the marked signature compared to the 
initial signature 

The correspondence rate between the initial 
signature and the detected signature remains to be 

15 defined. To do so, we perform the correlation 

A 

calculation between the two signatures S and S using 
the following operation: 

««*>- J***:™ 

A 

The two signatures correspond if cc(S,S)>T , 
20 where T is the decision threshold of the 
correspondence of the signatures. 

This correlation calculation allows to 
define the resemblance rate between the original 

A 

signature S and the detected signature S . In this 

25 particular embodiment of the invention, we set the 

threshold T to 0.7, so that we consider that if 

cc(S,S)> 0,7 , the detected signature corresponds to the 
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original signature . 

As previously indicated, the figures 8a and 
8b show an original image, and the corresponding 
marked image, according to the parameters of the 
5 particular embodiment of the invention shown above. 

Figures 9a to 9c illustrate, in the form of 
graphs, the results of tests of robustness of the 
mark inserted into the "home" image in figures 8a and 
8b. 

10 The graphs in figures 9a to 9c respectively 

correspond to the tests of the mark under attack from 
JPEG compression, from noise addition and from median 
filtration, and represent the mean of 20 measurements 
(use of 20 different keys K) for each attack. As can 

15 be seen, the results are stable, and the robustness 
of the watermarking technique of the invention 
against these attacks is very good. 

Although the present invention has been 
described with reference to preferred embodiments, 

20 workers skilled in the art will recognize that 
changes may be made in form and detail without 
departing from the spirit and scope of the invention. 
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APPENDIX 1: REMINDER ON THE THEORY OF WAVELETS 

1. Introduction 

The theory of wavelets allows to 
approximate a function (curve, surface, etc.) with 
5 different resolutions. Thus, this theory allows to 
describe a function in the form of a rough 
approximation and of a series of details allowing to 
perfectly reconstruct the original function. 

Such a multi-resolution representation 

10 (Mallat S., "A Theory for Multiresolution Signal 
Decomposition: the Wavelet Representation" , IEEE 
Transactions on Pattern Analysis and Machine 
Intelligence, Vol. 11, No. 7, July 1989, pp. 674- 
693.) of a function thus allows to interpret in a 

15 hierarchic manner the information held in the 
function. To do so, this information is reorganized 
and a set of details appear with different 
resolutions. Starting from a sequence of increasing 
resolutions (^j) JeZ / the details of a function with a 

20 resolution r. are defined as the difference of 
information between its approximation with the 
resolution r. and its approximation with the 

resolution r J+l . 

2. Annotations 

25 Before presenting the foundation of the 

multi-resolution analysis in fuller detail, we will 
present the annotation used in this section. 

- The sets of integral and real numbers are 
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respectively marked Z and R. 

- L 2 (R) designates the vectorial space of the 
measurable and integrable one-dimensional functions 
f(x) . 

5 - For f(x)eL 2 (R) and g(x)eL 2 (R), the scalar 

product of f (x) and g(x) is defined by: 

(f(x),g(x)}= fj(u)g(u)du . 

- For f(x)eL 2 (R) and g(x)e L 2 (R) , the convolution 
of f (x) and g(x) is defined by: 

10 f*g(x)=£f(u)g(x-u)du. 

- L 2 (R 2 ) designates the vectorial space of the 
functions f(x,y) of two measurable and integrable 
variables . 

- For f(x,y)eL 2 (R 2 ) and g(x,y) e L 2 (R 2 ) , the scalar 
15 product of f(x,y) and g(x,y) is defined by: 

(f( x > y)> g( x> y)) = £ %*J( u > v )g( u > v)dudv . 

3. Properties of the multi-resolution analysis 

In this section, we present in an intuitive 
manner the properties desired by the operator 

20 allowing to perform a multi-resolution analysis of 
the function. These properties are drawn from work by 
Mallat S., in yy A Theory for Multiresolution Signal 
Decomposition: the Wavelet Representation", IEEE 
Transactions on Pattern Analysis and Machine 

25 Intelligence, Vol. 11, No. 7, July 1989, pp. 674-693. 

Let A j be the operator which approximates a 
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function f(x)sL 2 (R) with the resolution 2 J (j>0) (i.e. 

f (x) is defined by 2 J samples). 

The expected properties of A 2i are the 

following : 

5 1. A 2 j is a linear operator. If A 2j f(x) represents 

the approximation of f (x) with the resolution 2 J , 
then A 2 jf(x) must not be modified when a new 

approximation is made with the resolution 2 J . This 
principle is written AjoAj=A 2J and shows that the 

10 operator A 2J is a projection operator in the vectorial 

space V 2 jCzL 2 (R). This vectorial space can be 
interpreted as the set of all the possible 
approximations with the resolution 2 J of the function 
of L 2 (R). 

15 2. Among all the possible approximations of f (x) 

with the resolution 2 J , A 2J f(x) is the most similar to 

f (x) . The operator A^ is therefore an orthogonal 

projection onto V 2J . 

3. The approximation of a function with the 

20 resolution 2 y+1 contains all the information necessary 
to calculate the same function with the lower 

resolution 2 J . This casualty property induces the 
following relation: 
VjeZ,V 2j czV 2j+l . 

25 4. The approximation operation is the same with 
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all the resolutions. The spaces of the approximated 
functions can be derived one from another by a change 
in scale corresponding to the difference in 
resolution: 
5 V/ e Z, f(x) € V 2J « f(2x) e . 

5. When we calculate an approximation of f (x) 

with the resolution 2 J , a part of the information 
contained in f (x) is lost. However, when the 
resolution tends towards the infinitive, the 
10 approximated function must converge towards the 
original function f(x). In like manner, when the 
resolution tends towards zero, the approximated 
function contains less information and must converge 
towards zero. 

15 Every vectorial space (Vj) JeZ which satisfies 

all of these properties is called approximation 
multi-resolution of L 2 (R) . 

4. Multi-resolution analysis of a one-dimensional 
function 

20 4.1 Search for a base of V t 

As is indicated above, the approximation 
operator A s is an orthogonal projection onto the 

vectorial space V 2j . In order to numerically 
characterize this operator, an orthogonal base of V v 
25 must be found. being a vectorial space containing 

the approximations of the functions of L 2 (R) with the 
resolution 2 J , every function f(x)eVj can be seen as 
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a vector with 2 j components. Therefore, 2 j basic 
functions must be found. 

One of the principal theorems of the theory 

of wavelets stipulates that there is a single 

5 function 0(x)eL 2 (R) , called the scale function, from 
which we can define 2 j basic functions Of(x) of Vj 
via dilatation and translation of O(x) : 

Of(x ) = 0(2 J x - i), i = Or • -,2 y -1. 
Approximate a function f(x)eL 2 (R) with the 
10 resolution 2 J therefore returns to project in an 
orthogonal manner f (x) onto the 2 j basic functions 
Of (x) . This operation consists in calculating the 

scalar product of f (x) with each of the 2 J basic 
functions 0\(x) : 

k=2 J -\ 



15 



A 2J f(x)= Z(f(u),Oi(u))oi(x) 

= k 'f\f(u),<D(2 J u-kffi>(2 J u- k). 



k=Q 

We can demonstrate that Ajf(x) can reduce 
the convolution of f (x) with the low-pass filter O(x) , 
evaluated at point k : 

A 2j f = (f(u)*0(-2 J u))(k) f keZ. 
20 As O(jc) is a low-pass filter, A s f can be 

interpreted as a low-pass filtering followed by an 
even under-sampling . 
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4.2 Construction of the multi-resolution analysis 

In practice, the approximated functions f 
(signal, image, etc.) are discrete. Supposing that 
the original function f (x) is defined on samples 
5 n = 2 k (keZ) . The maximum resolution of f (x) is then n. 

Let A n f be the discrete approximation of 
f (x) with the resolution n. The casualty property 
(see §3 above) claims that A,-f can be calculated from 

A n f for every j<k. 
10 Indeed, by calculating the projection of 

the 2 J basic functions ®l(x) of V 2j onto F 2 , +1 , it can 

be demonstrated that A jf can be obtained by 

convoluting A^f with the low-pass filter 

corresponding to the scale function and by under- 
15 sampling the result of a factor 2: 

A 2J f(u) = 2 'j?h(k - 2u)A 2j+l f(k),0 <u<2 J -\ 

with h(n) = (®(2u),<b(u-n))yneZ . 

4.3 The detail function 

As is mentioned in property (5) of §3, the 
20 operation which consists in approximating a function 

f (x) with a resolution 2 J from an approximation with 

the resolution 2 ;+1 generates a loss of information. 
This loss of information is contained in a function 

called a detail function with the resolution 2 j and 
25 designated D f . It is to be noted that the knowledge 
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of D 2 jf and A^f allows to perfectly reconstruct the 
approximated function A 2j+l f . 

The detail function with the resolution 2 J 
is obtained by projecting in an orthogonal manner the 
5 original function f (x) onto the orthogonal complement 
of V 2 j in V . Let W j be this vectorial space. 

To numerically calculate this projection, 
an orthogonal base of Wj must be found, i.e. 2 s basic 

functions. Another major theorem of the theory of 
10 wavelets stipulates that from a scale function O(x) , 

it is possible to define 2 s basic functions of W 2i . 

These basic functions ^/(x) are obtained by 
dilatation and translation of a function *¥(x) called 
wavelet function: 
1 5 Vf(x ) = *¥(2 J x - i), i = 0, - • -,2 y - 1. 

In the same manner as for the construction 
of the approximation A 2j f , we can demonstrate that 

D jf can be obtained by convoluting the original 

function f (x) with the high-pass filter *¥(x) followed 

20 by an under-sampling of a factor 2 J : 

D 2i f = (f(u)*V(-2 J u))(k),keZ. 

4 . 4 Extension of the multi-resolution analysis of 
two-dimensional functions 

In this section, we present the method of 
25 extending the multi-resolution analysis by wavelets 

to the L 2 (R 2 ) functions of as the images. To do so, 
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we use the same theorems as those previously used. 
Thus, if we designate V v the vectorial space of the 

approximations of L 2 (R 2 ) with the resolution 2 y f we 
can demonstrate that an orthonormal base of can be 
found by dilating and translating a scale function 
0(x t y)eL 2 (R 2 ): 

®l(x,y) = 0(2 J x - i,2 J y - j), (ij) e Z 2 . 

In the particular case of separable 
approximations of L 2 (R 2 ), we have ®(*>y) = ®(*my) 
where <S>(x) is a scale function of L 2 (R) . In this case, 
the multi-resolution analysis of a function of L 2 (R 2 ) 
is carried out by processing in a sequential and 
separable manner each of the dimensions x and y. 

As in the one-dimensional case, the detail 

function with the resolution 2 y i s obtained by an 
orthogonal projection of f(x,y) onto the complement 
of V 2 j in V 2j+l , designated W j . In the two-dimensional 

case, we can demonstrate that if we designate *F(x) 
the wavelet function associated with the scale 
function O(x) , then the three functions defined by: 

V l (x,y) = Q(x/¥(y) 
V 2 (x,y) = V(x)®(y) 
^(x,y) = V(x)V(y) 

are wavelet functions of L 2 (R 2 ) . By dilating and 
translating these three wavelet functions, we obtain 
an orthonormal base of : 
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(x,y) = ^(V x-k.Vy-l) 
x ¥ 2 (x,y) = V®(2 J x-k,2 J y-l) 
*¥ i (x,y) = W(Vx-k,2 i y-l). 

The projection of f(x,y) onto these three 
functions of the base W , give three detail functions: 

D l 2j f = f(x,y)*<!> J (-x)V j (-y) 
D 2 2i f = f(x,y)*V J (-x)<S> j (-y) 
D 2j f = f(x,y)*V J (-x)Vj(-y) 



